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We have studied the tunneling conductance in two dimensional electron gas / insulator / super- 
conductor junctions in the presence of Rashba spin-orbit coupling (RSOC). It is found that for low 
■ insulating barrier the tunneling conductance is suppressed by the RSOC while for high insulating 

barrier it is almost independent of the RSOC. We also find the reentrant behavior of the con- 
ductance at zero voltage as a function of RSOC for intermediate insulating barrier strength. The 
• results are essentially different from those predicted in ferromagnet / superconductor junctions. The 

present derivation of the conductance is applicable to arbitrary velocity operator with off-diagonal 
components. 
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I. INTRODUCTION 



In normal metal / supercunductor (N/S) junctions Andreev reflection (AR) 1 is one of the most important process 
for low energy transport. The AR is a process that an electron with up spin injected from N at the energy below 
^ ■ the energy gap A is converted into a reflected hole with up spin. To describe the charge transport in N/S junctions 
Blonder, Tinkham and Klapwijk (BTK) proposed the formula for the calculation of the tunneling conductance 2 . 
A gap like structure and the douling of tunneling conductance appear in the voltage dependence due to the AR. 
This method has been extended to the ferromagnet / superconductor (F/S) junctions and used to estimate the spin 
\ polarization of the F layer experimentallji^i^. In F/S junctions, AR is suppressed because the retro-reflectivity is 
broken by the exchange field in the F layer—. As a result, the conductance of the junctions is suppressed-. Spin 
dependent transport in F/S junctions is an important subject in the field of spintronics which aims to fabricate novel 
i devices manipulating electron's spin. 

Spintronics has recently received much attention because of its potential impact on electric devices and quantum 
t— ( \ computing 8 . Among recent works, many efforts have been devoted to study the effect of spin-orbit coupling on 
J> ■ transport properties of two dimensional electron gas (2DEG)2*i2*ii*i2*i2*i^. The pioneering work by Datta and Das 
suggested the way to control the precession of the spins of electrons by the Rashba spin-orbit coupling (RSOC)i£ in 
F/2DEG/F junctions^. This spin-orbit coupling depends on the applied electric field and can be tuned by a gate 
' voltage. On the other hand spin dependent transport based only on spin-orbit coupling without ferromagnet, e.g., 
, spin Hall effect is also a hot topi o 17 i 18 ' 19 . 

As in the case of exchange field in F/S junctions, RSOC may affect the tunneling conductance in 2DEG/S junc- 
tions because RSOC mixes spin-up and spin-down states. The RSOC induces an energy splitting which lifts the 
spin degeneracy, but the energy splitting doesn't break the time reversal symmetry unlike an exchange splitting in 
ferromagnet (see Fig. QJ. Therefore transport properties in 2DEG/S junctions may be qualitatively different from 
those in F/S junctions. However, in 2DEG/S junctions the effect of RSOC on transport phenomena is not studied 
well. It is desirable to make a formalism incorporating the effect of the RSOC in these junctions. For this purpose a 
BTK-like formula may be accessible. However the derivation of the conductance by BTK cannot be directly extended 
to that in 2DEG/S junctions because velocity operator has off-diagonal components by RSOC. 
£j , In this paper we present a general method to derive a conductance in superconducting junctions which is applicable 
IV ' to arbitrary velocity operator with off-diagonal components. Applying it, we calculate the tunneling conductance in 
. £h ! 2DEG/S junctions, compare it with that in F/S junctions and clarify how RSOC affects the AR and normal reflection 
probabilities. The obtained results can be useful for the design of mesoscopic 2DEG/S junctions and for a better 
understanding of related experiments. 

The organization of this paper is as follows. In section II, we will provide the detailed derivation of the expression 
for the conductance. In section III, the results of calculations are presented for various types of junctions. In section 
IV, the summary of the obtained results is given. In the present paper we confine ourselves to zero temperature. 



II. FORMULATION 



We consider a ballistic 2DEG / S junctions where the 2DEG/S interface is located at x = (along the y-axis), and 
has an infinitely narrow insulating barrier described by the delta function U{x) — U5{x). 
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FIG. 1: (color online) Schematic illustration of Rashba and Zeeman splitting. 
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FIG. 2: (color online) Schematic illustration of scattering processes. 



The effective Hamiltonian with RSOC is given by 
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with k± = k x ± iky, the energy gap A, = £— (k 2 — tip), Fermi wave number kp, Rashba coupling constant A, and 
step function 9{x). 

Velocity operator in the x-direction is given by2£ 
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As shown in Fig. [2 the wave function ip(x) for x < is represented using eigenfunctions of the Hamiltonian: 



ip(x < 0) = e lk y y 
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for an injection wave with wave number fc^) where k\ = — + ^/ (^r) 2 + ^fj 

^i(2)± — fci(2) e± * ei<2> ■ a i(2) and ^1(2) are AR coefficients. Cx(2) and dx(2) are normal reflection coefficients. #1(2) is an 
angle of the wave with wave number ku2) with respect to the interface normal. 
Similarly for x > ift(x) is given by 



ip{x > 0) = e ifc ^ 
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where -E is quasiparticle energy and 9 is an angle of the wave with wave number kp with respect to the interface 
normal. e-i<z), /i(2)>ffi(2) and /ii(2) are transmission coefficients. Note that since the translational symmetry holds for 
the ^-direction, the momenta parallel to the interface are conserved: k y = kF sin# = k\ sin#i = k2 sin 62- 
The wave function follows the boundary conditions 2 ^: 
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Now we will derive a formula for the tunneling conductance. Before giving detailed calculation, we present an 
essential idea for the derivation of the conductance: First we calculate expectations of current for the complete sets of 
the eigenfunctions. Next we sum the expectations multiplied by corresponding distribution functions. Then we can 
get the net current of the junctions. Detailed derivation is given in the Appendix. 

Finally we obtain the dimensionless conductance represented in the form: 
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JVi and N2 are density of states normalized by those with A = for wave number k\ and hi respectively. A12 and A21 

are defined in the Appendix. The critical angle 9c is defined as cos 9c = J 

<jn is given by the conductance for normal states, i.e., as for A = 0. We define normalized conductance as 
a T = as/a N and parameters as f3 — jrgj^- and Z = fi|^-. For example, in InGaAs heterostructures, (3 is estimated 

as (3 ~ 0.2m±£l Here we choose the same effective mass in 2DEG and S. In most cases the effective mass in 2DEG is 
much smaller than that in S. However we can show that this effect is equivalent to that by the increase of Z. Thus we 
neglect the difference of the effective masses in the present paper. 



III. RESULTS 



First we study the normalized tunneling conduntace ot as a function of bias voltage V in Fig. [3] For Z = 10 where 
the AR probability is low, ctt is almost zero within the enregy gap and independent of j3. In contrast, for Z = 1, ctt 
is slightly enhanced with the increase of (3 around zero voltage. For Z = where the AR probability is very high, <tt 
becomes two for /3 = within the energy gap. It is reduced by the increase (3 within the energy gap. 




FIG. 3: (color online) Normalized tunneling conductance with Z = 10 in (a), Z — 1 in (b), and Z = in (c). 

Next we study the difference between the effect of the Rashba splitting and that of Zeeman splitting. We have 
calculated conductance in F/S junctions following Ref. 7 . We plot the tunneling condutance for superconducting states 
as at zero voltage for 2DEG/S junctions in (a)-(c) and for F/S junctions in (d)-(f) of Fig. |2]with Z = 10 in (a) and 
(d), Z = 1 in (b) and (e), and Z = in (c) and (f). In (a)-(c) we show the dependence of as, normalized by ctat for 
(3 = 0, on (3 for various Z . For Z = 10 it has an exponential dependence on [3 but its magnitude is very small while 
it has a reentrant behavior as a function of (3 for Z = 1. For Z = it decreases linearly as a function of (3. On the 
other hand, in F/S junctions, the dependence of as on U, normalized by Fermi energy Ep, is qualitatively different. 
We plot a s normalized by <7jv at U = 0. The exchange field suppresses as independently of Z as shown in (d)-(f). 
This is because the AR probability is reduced by the exchange field. Therefore the effect of the Rashba splitting on 
conductance is essentially different from that of Zeeman splitting on conductance. This can be explained as follows. 
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The Zeeman splitting gives unbalance of populations of up and down spin electrons. Thus it suppresses the AR where 
pairs of spin-up and spin-down electrons are transmitted to S. On the other hand, the Rashba splitting never causes 
such an unbalance. Thus it cannot suppresses the AR, which results in various (3 dependence of the conductance. 




P U IE F 



FIG. 4: Tunneling condutance for superconducting states at zero voltage as a function of RSOC in 2DEG/S junctions (left 
panels) and the exchange field in F/S junctions (right panels) with Z — 10 in (a) and (d), Z — 1 in (b) and (e), and Z = in 
(c) and (f). Here /? = 

In order to explain the line shapes of the conductances, we will check the angular averaged normal reflection and AR 
probabilities as a function of voltage. For large Z, AR probabilities are small and normal reflection probabilities reflect 
the dependence of the densities of states on f3: N% is a decreasing function of (3, while N2 is an increasing function 
of j3. Therefore C\ and C2 are reduced and D\ and D2 are enhanced with the increase of (3. Figure [5] shows the 
probabilities for Z = 10. AR probabilities {A\, A2, B\ and B2) are slightly enhanced around eV = A and have similar 
structures with the increase of [3 while normal reflection probabilities D\ and D2 increase with the increase of (3. On 
the other hand normal reflection probabilities C\ and C2 arc reduced with the increase of (3. In other words, reflected 
waves with wave number fci(fc2) are suppressed (enhanced) by RSOC. The enhancement and the suppression compete 
with each other. Thus the conductance is almost independent of RSOC. For small Z, normal reflection probabilities 
are small. AR probabilities A\ and B2 are reduced as increasing (3. This stems from the mismatch of Fermi surfaces 
between 2DEG and S by the increase of (3. In fact, for Z = (see Fig. |H} normal reflection probabilities C2 and D\, 
and AR probabilities B\ and B2 are slightly enhanced with the increase of (3. Normal reflection probabilities C\ and 
D2 increase by the increase of (3. On the other hand AR probabilities A\ and B2 within the energy gap are reduced 
with the increase of (3. This means that only eigenfunctions with the same wave number as the injection wave are 
affected by RSOC. From Fig. we can understand the suppression of the tunneling conductance by RSOC (see Eq. 
(7))- 
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FIG. 6: (color online) The angular averaged Andreev and normal reflection probabilities for Z = 0. A\, A2, B\ and B2 are AR 
probabilities. Ci, C2, D\ and D2 are normal reflection probabilities. 
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IV. CONCLUSIONS 

In the present paper we have studied the tunneling conductance in two dimensional electron gas / insulator / 
superconductor junctions with RSOC. We have extended the BTK formula and calculated the tunneling conductance. 
It is found that for low insulating barrier the tunneling conductance is suppressed by the RSOC while for high 
insulating barrier the tunneling conductance is almost independent of it. We also found a reentrant behavior of the 
conductance at zero voltage as a function of RSOC for intermediate insulating barrier strength. This phenomena are 
essentially different from those found in F/S junctions where the tunneling conductance is suppressed by exchange 
field, being independent of the barrier strength. The present derivation of the conductance is applicable to arbitrary 
velocity operator with off-diagonal components. 

The results give the possibility to control the AR probability by a gate voltage. We believe that the obtained results 
are useful for the design of mesoscopic 2DEG/S junctions and for a better understanding of related experiments. 

In this paper we focus on ballistic 2DEG/S junctions. In diffusive 2DEG/S junctions, proximity effect plays an 
important role. The RSOC breaks the inversion symmetry and hence mixes the parity. As a result, "triplet" pairing 
may be induced in the 2DEG region^ as predicted in diffusive F/S junctions 2 —. The study in this direction is now in 
progress. 
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Nanoscience Project, the Ministry of Education, Culture, Sports, Science and Technology, Japan, the Core Research 
for Evolutional Science and Technology (CREST) of the Japan Science and Technology Corporation (JST) and a 
Grant-in- Aid for the 21st Century COE "Frontiers of Computational Science" . The computational aspect of this 
work has been performed at the Research Center for Computational Science, Okazaki National Research Institutes and 
the facilities of the Supercomputer Center, Institute for Solid State Physics, University of Tokyo and the Computer 
Center. 



Appendix 

Here we give a detailed derivation of the conductance which is applicable to arbitrary velocity operator with 
off-diagonal components. For an electron injection and a hole injection from 2DEG (represented by ip e and iph, 
respectively), the resulting currents j e and jh in the 2DEG region are given by 

i e = ReOKrg^e) oc (l + A he + B he - C ee - D ee ) (9) 

Jh = Re(^T 3 ^) « (1 + A eh + B eh - C hh - D hh ) . (10) 

Similary, for an electron and a hole injection from S (represented by ipi e and ipfh, respectively), the corresponding 
currents j' e and j' h in the 2DEG region reads 

j' e = Re(V»?»«7sO cx (F ee + G ee - H he - J he ) (11) 



j' h = Re«W 3 ^) oc (F hh + G hh - H eh - J eh ) . (12) 

Here A he and B he , and, C ee and D ee denote AR and normal reflection probabilities with electron injection respectively. 
F ee and G ee are normal transmission probabilities with electron injection. H he and J he are transmission probabilities 
with the injection of an electron converted into a hole at the interface. Other notations are defined in a similar 
way. Note that there are four independent eigenfunctions: two kinds of electron- like quasiparticles and two kinds of 
hole-like quasiparticles. The total current reads 

/oo 
( Je f(E - eV) - Jh f(E + eV) + j'J(E) - j' h f(E))dE 
-OO 



/OO 
(1 + A'l e + B\ e - Cf - Df) (f(E - eV) - f(E + eV))dE. 
-oo 



(13) 
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Then the differential conductance at zero temperature has the form 

dl 



dV 



oc (l + A he + B he - C ee - D ee ) (14) 



with Fermi distribution function f(E) and bias voltage V. Here we assume the particle-hole symmetry which results 
in the relations X ee = X hh and X he = X eh (X = A,B,C,D,F,G,H,J). The original BTK method^ cannot treat 
velocity operator with off-diagonal components. However, the derivation given here is applicable to arbitrary velocity 
operator with off-diagonal components. 

Let us apply the above procedure to our model. For an injection wave with wave number fci, the current reads 

il = ^ (l + Q fci cosfc + | ai | 2 + \b x \ 2 ^A 2 i - | Cl | 2 - \d x \ 2 ^A 2 i) . (15) 
For an injection wave with wave number k2, the current is 



2m 

Here we define A12 and A21 as 



= ^ f 1 + y J fc 2 cos0 2 f 1 + |a 2 | 2 ^A 12 + |6 2 | 2 - |c 2 | 2 ^A 12 1 ) . ( Ki) 



fci cos Ox k 2 cos8 2 , . 

A12 — 7- A21 = t-. (17) 

fc2 COS #2 fci COS #1 
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